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ABSTRACT 

We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy 
medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and 
have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent 
three dimensional radiative transfer code. We provide an AGN library of spectral energy distributions (SEDs). Its purpose is to quickly 
obtain estimates of the basic parameters of the AGN, such as the intrinsic luminosity of the central source, the viewing angle, the inner 
radius, the volume Ailing factor and optical depth of the clouds, and the optical depth of the disk midplane, and to predict the flux at 
yet unobserved wavelengths. The procedure is simple and consists of finding an element in the library that matches the observations. 
We discuss the general properties of the models and in particular the 10 pm silicate band. The AGN library accounts well for the 
observed scatter of the feature strengths and wavelengths of the peak emission. AGN extinction curves are discussed and we And that 
there is no direct one-to-one link between the observed extinction and the wavelength dependence of the dust cross sections. We show 
that objects of the library cover the observed range of mid IR colors of known AGN. The validity of the approach is demonstrated 
by matching the SEDs of a number of representative objects: Four Seyferts and two quasars for which we present new Herschel 
photometry, two radio galaxies, and one hyperluminous infrared galaxy. Strikingly, for the five luminous objects we And pure AGN 
models fit the SED without a need to postulate starburst activity. 
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1. Introduction 


The suggestion that the 2.2-22pm emission of Seyfert galax¬ 
ies and quasars may be due to d ust ([Low & Kleinmann|p^68] ) 
was first explored analytically by [Rees et ( |1969| ) who were 
able to reproduce the slope of the near- to mid IR spectrum. Ra¬ 
diative transfer models of spherically symmetric distributions of 
dust in AGN were first developed by |Rowan-Robins on & Craw¬ 


ford] (|1989|), [Lawrence et al.| ( ri991| ), and |Rowan-Robinson & Ef- 
stathiou| ( |1993| ). 


Trying to account for the infrared emission of type 1 and 
type 2 AGNs with models that assume a disc-like or toroidal 
distribution of dust was recognized e arly on as an im portant 
test of the unified model for AGN ( |AntonuccI||1993|). Such 
models have been developed from the early 1990s ([Efstathiou 


& Rowan-Robinson||1991t |Pier & Krolik||1992 1993t [Rowan 


Robinson et al.||1993r Granato & Danese||1994t [Efstathiou & 

Rowan-Robinsonl 1 9^ Efstathiou et al.|1995 1. Most of the dis¬ 
cussion on these models concerned the interpretation of the ob¬ 
served behaviour of the silicate feature at 9.7 pm which was ab¬ 
sent in emission in type 1 objects ( [Roche et al.|1991| ) and chal¬ 
lenged the whole idea of trying to explain the infrared emis¬ 
sion of AGN with dust models and the unified model itself. It 


is now recognized that smooth and clumpy torus models (Eel- 
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* Herschel is an space observatory wi tn scien ce instruments pro¬ 

vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 


tre et 

al.|2012 ) can reproduce the whole range of features which 

are o 

^served using Spitzer spectroscopy (Siebenmorgen et al. 

2005 

[Hao et al. [2005[[Spoon et al.[2007[ [Schweitzer et al.[2008[ 

Menc 

oza-Castrejon et al.[[2015l[Hatziminaoglou et al.[[2015|) as 

well as high spatial resolution ground based data (Honig et al. 

2010l[Alonso-Herrero et al.|2011l[Gonzalez-Martm et al.|2013l 

Esquej et al.|2014l Ramos Almeida et al.|2014a[). Another chal- 


lenge for these models is trying to explain the spatial extent of 
the mid IR emission which also looked puzzling initially. A sig¬ 
nificant part of the mid IR emission appears to come from the 


ionization cones (Braatz et al.[ 1993 [[Cameron et al.[ 1993 [[Honig 

etal.2012 2013 

), which suggests the presence of dust in the ion- 

ization cones eit 

ler in the form of discrete clouds or as outflows 

([Efstathiou et al. 

1995[[Elitzur & Shlosman[2006[). The data sug- 

gest that what is 

needed might be a combination of smooth and 


clumpy structures in AGN which may extend into the narrow 
line region, see [Netzerj ( |2015| ) for a recent review. The present 
picture is also supported by recent hydrodynamical simulations 
which predict a filamenta ry torus and a turbulent disc component 
( jSchartmann et al.|2Q14| ). 


Large ground-based telescopes offer an opportunity to much 
better isolate the contribution of the host galaxy and circumnu- 
clear star formation from the AGN when compared to observa¬ 
tions from spacejAlonso^HOTero^£F^|20^ [Ramos Almeida 


et al.[[20lT| [Asmus et al.[[2Q14| ). Spitzer spectra provide a res 
olution of ~ 3” whereas high resolution MIR spectra from 
the ground are an order of magnitude sharper, and sensitive 
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enough to observe local Seyferts at luminosities < 10"^^ erg/s. 
Such high resolution studies have been undertaken for exam- 
pie in NGC 1068, NGC 1365, NGC 3281, and Circinus (iMa^ 


son et al.|2012| ^che et al.|2006t|Alonso-Herrero et al.|2012t 


Sales et al. mm Esquej et al.||2014|) A sample of 19 Seyferts 

are presented by [Honig et al. ( 201()|) and another 22 Seyferts 
are st udied by |Gonzalez-Martfn et al.| ( [20T3| ). At ~ 20 pc reso¬ 
lution |RuscherDutoetar[^20T4]rd^cted the silicate feature at 
9.7 yum in absorption for the Seyfert 2 NGC 1386, and in emis¬ 
sion for t he Seyfert 1 NGC 7213 with a p eak wavelength beyond 
10.5 yum. [Ramos Almeida et al.| ( |2014a| ) compared the nuclear 
MIR spectrum with the Spitzer spectrum and find that the AGN 
in Mrkl066 dominates the continuum emission at d < 15yum 
on scales of ~ 60 pc (90 - 100%) whereas this contribution de¬ 
creases to below 50% when the emission of the central ~ 830 pc 
is considered. 

IR interferometry offers unprecedented angular resolution 
that a l lows resolving the inner regi ons of dusty AGN ( jJaffe et ah 
|2004[ [Tristram et al.||2007] [2009| ). It has been shown that the 
inner torus radii derived from visibility curves as expected scale 
approximately with the square root of the A GN luminosity ( [Sug- 
[anuma et aH2006[ [Kishimoto et al.|2()TT] ), except for some de¬ 
viations that are observed for fainter sources ( [Burtscher et al.| 


2013| ). A sam ple of nearby Seyfe rt galaxies observed with IR in¬ 
terferometers ( [Honig et al.|2013[ ) reveal a diversity of obscuring 
structures. A significant part of the mid IR emission of certain 
sources appears to come from the the polar direction. These ob¬ 
servations are difficult to reconcile with torus models alone and 
challenges the justification of using clumpy models to fit that 
emission ( [Honig et al.|2013| ). Polar dust may be physically dis¬ 
connected from the torus and may arise from a radiatively driven 
dusty wind blown off the inner region of the torus. 

ALMA observations have recently been reported for the 
Seyferts NGC 1566 ( [Combes et al.[2014]), N GC 1068 ( [Garcia-[ 
[Burillo et al.[[2014[ ) and NGC 34 ( [Xu et al.[[2Q14[ ). Typically at 
25-50pc resolution an unresolved core and a circumnuclear disk 
extending about 200pc in radius is detected. 

Solutions of the radiative transfer in clumpy AGN tori em¬ 
ploy methods that have diff erent levels of sophisti cation. A first 
formalism was presented by [Nenkova et al.[ ( [2002[ ) utilizing aid 
radiative transfer code and computed SED of clouds by assum¬ 
ing a slab geometry. The slabs are heated from different angles 
and then averages are taken so that they mimic the emission of 
spheres. These spheres are arranged so that a torus-like geome¬ 
try is sim ulated. [Honig et al.[([2006[) pr e-computed these clouds 
in 2d and [Honig & Kishimoto[ ( [201Q| with an upgrade of that 
model in 3d. In the AGN torus there are clouds that are not il¬ 
luminated directly by the central source because other optically 
thick clouds lie along the line of sight towards the primary heat¬ 
ing source. Such clouds are shadowed and indirectly heated. In 
the above mentioned models it is assumed that this indirect heat¬ 
ing is coming from clouds in the immediate surrounding which 
are directly heated so that averages can be build. Clumps, which 
are directly heated by photons reflected on a surface of an opti¬ 
cally thick disk were not treated. Clumpy torus models are suc¬ 
cessful in explaining the strengths of the silicate emission and 
absorption band. However, they do not account for ver y deep sil- 
icate absorptions and under-predict the NIR emission ( [Nenkova 
[et al.|2010 [Spoon et al.|2007[). Self-consistent 3d clu mpy torus 
^ y [Schart 


models are presented by|Schartmann et al.[([2008| ) and [Stalevski[ 


|et al.[ ( [MT2a| ). [Mendoza-Castrej6n et al.[ ( [2Q15[ ) find that Seyfert 


galaxies have clumpy dust distributions and/or a disk compo¬ 
nent. The formation of a clumpy and a continuous dust compo¬ 
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nent is also pred icted by hydrodynamical simulations ( [Schart- [ 
[mann et al.|2Q14[ ). 

In this paper we describe our SED library of 2-phase 
AGN torus models that has been computed with the fully self- 
consistent 3d Monte Ca rlo radiative transfer code of [Heymami] 
[& Siebenmorgen[ ( [20 1 2 ) . We apply beside the clumps a density 
distribution of a disk component that extends into the ionization 
cones. We will show below that the emission from such diffuse 
dust in the ionization cones can dominate the total emission from 
an edge-on view. This can explain why some Seyferts show ex¬ 
tended emission coincident with the ionization cones at lOpm 
despite the fact that there is an optically thick torus which lies 
in a plane perpendicular to the ionization cones. The emission 
coming from the ionization cones is also much more compact 
than the torus in agreement with interfer ometric observations 
( [Burtscher et aL|2()l3[ [Honig et aL|20T3| l. We discuss the gen¬ 
eral properties of the models (Sect. 3) and apply the SED library 
to fit a number of representative objects (Sect. 4). We summarize 
our findings in the conclusions (Sect.5). 



Fig. 1. Schematic view of the AGN structure assumed in this work. 
Dust extends from the dust evaporation radius at pc/10 < rin < pc/2 up 
to distance of rout = 170 X Pn from the black-hole (BH). Note that dust 
from the smooth disk emerges also as a diffuse component into the polar 
region. Dust clouds are pc sized with cube length of J = 3 x rin. The hot, 
warm and cold side of the clouds is indicated and for one clump their 
shadows. The accretion disk in the immediate black-hole environment 
and the jet is not treated. 


2. AGN Model 

In this section we describe the AGN model library of SEDs. We 
outline the radiation transfer code, the assumed dust structure, 
and the library set-up. We choose a condensed description of the 
AGN model with all parameters kept constant except a minimum 
of five key parameters that are varied to compute the library. We 
visualize the dust density distribution and AGN images at differ¬ 
ent wavelengths, and show the impact of our free parameters on 
the SED. 

Most models for the torus in AGN consider dust that is rep¬ 
resentative for the diffuse interstellar medium and the fact that 
there are huge differences in the density and the radiation envi¬ 
ronment of both media is ignored. Our library is computed con¬ 
sidering fluffy grains. The rational for this and the dust proper¬ 
ties assumed are discussed in Sect. [2.6[ Einally, we discuss the 
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impact of scattering in spatially unresolved observations of the 
extinction curve in AGN. 


2.1. Radiative Transfer 


For our calculations we apply the Monte Carlo (MC) continuum 
radiative transfer code of [Heymann & Siebenmorgen] ( |2012| ). 
The code allows us to follow the light path of photons that 
are emitted from a heating source through an arbitrary three 
dimensional (3d) dust conhguration. Photons are binned into 
packets of equal energy, those photons of one packet have the 
same frequency. Dust absorption, scattering and emission events 
are computed by chance. The process depends on the optical 
depth and the absorption and scattering cross section of the dust. 
We utilize the MC algorithm developed by [Bjorkman & Wood] 
( [200 1| ), which is iterative-free as it employs an immediate up¬ 
date of the dust temperatures. For very optically thin regions we 
use the method by jLucyj ( |1999| ) which reduces uncertainties in 
the dust temperatures. Anisotropic scattering is treated assuming 


the phase function p(cos(6)) as given by (Henyey & Greenstein 

[T9^ 


p(COS(e)) : 




4n (I -2g cos(6))^i^ 


( 1 ) 


where 6 is the scattering angle. The anisotropy factor g of the 
grains is provided by the dust model and whenever g = 0 we 
consider isotropic scattering. 


The model space is set up as cubes in an orthogonal Cartesian 
grid. Each cube can be divided into an arbitrary number of sub¬ 
cubes with constant density. An illustration of the grid is given in 
Fig.l by Heymann & Siebenmorgen ( |2012| ). The subdivision of 
cubes allows a hner sampling whenever required. The advantage 
of such a grid is that it leads to a huge reduction of the computer 
memory for some applications. Examples of this are the interface 
of optically thick clumps that are embedded in an otherwise ho¬ 
mogeneous medium, regions close to the dust evaporation zone 
or conhgurations where the ratio of the outer and inner radius 
from the central source where dust exists is large rout/^n >100. 
The number of cells can be reduced further by a factor of eight 
by locating the source at origin and assuming mirror-symmetry 
with respect to the planes at x = 0, y = 0 and z = 0. Vectoring 
the code so that it can run on shared memory devices reduces 
the huge computing time of the MC procedure. The number of 
photon packets that can be followed in parallel on such machines 
scales roughly with the number of processing units available and 
reduces the required computing time accordingly. 


Ray-tracing in parallel projection from the detector of the ob¬ 
server through the model space, using temperatures and scatter¬ 
ing events for each cell of the MC run, allows the generation of 
high spatial resolution and high signal-to-noise images. In this 
work we derive the SED of the object by counting the photon 
packets that eventually escape the model space into a particular 
solid angle. This method gives within the photon noise the same 
results when compared to SEDs computed by ray tracing. The 
later method has the advantage that it provides images for each 
frequency. However, it is expensive in computing time. Eor one 
AGN model and SEDs computed in 9 viewing direction the ray 
tracer takes about half a day of computing time while the photon 
counting is completed in a few seconds. 


2.2. AGN Structure 


The super-massive central black hole with its accretion disk 
( Shakura & Sunyaevj |1976| ) is considered to be the pri¬ 
mary heating source of the circumnuclear environment. It 
emits anisotropic hard UV and optical photons ( |Netzer|[T987t 
[Kawaguchi & Mori|201 Ij ). The anisotropy of the primary source 
does not have a strong impact on the torus spectrum when com¬ 
pared to models with isotropic heating. The overall SED shape 
and the 10 pm Si-0 dust feature a re not seriously affected by the 
anisotropy ( [Stalevski et al.[2Ql^ . Therefore the emission of the 
primary heating source is taken to be isotropic. In optically thick 
applications the spectral shape of the heating source has virtu¬ 
ally no impact on the dust emission spectrum. We assume that 
the primary energy source takes a spectral shape approximated 
by [Rowan-Robinson| ( [1995[ ): 




for 

—oo 

< 

A < 0.01 

Oum) 


const 

for 

0.01 

< 

A <0.1 

(jum) 

1 oc < 

^-0.5 

for 

0.1 

< 

A < 1 

Oum) 



for 

1.0 

< 

A < -1-00 

Oum) 

The inner radius Tin 

of the dust torus is set close 


( 2 ) 


evaporation region, which scales with the total AGN luminosity 
as Lagn ^ ^ AGN this relation is conhrmed by re¬ 

verberation measurements and interferometry ( [Kishimoto et al.| 


2011). As long as this relation is respected, the SED becomes 


scale invariant, so that the SED shape does not vary for differ¬ 
ent AGN luminosities. We therefore keep the luminosity of the 
primary source constant at Ln = 10 ^^Lq. The sublimation tem¬ 
perature of the dust depends on details such as grain mineralogy, 
porosity and size. We assume that it is somewhere between 800 
- 1800 K and use Tin as a free-parameter. 

The required maximum size of the torus depends on the ob¬ 
served wavelength. However, choosing rout smaller than the size 
corresponding to the temperature, respectively wavelength, of 
interest results in an artihcial cut-off in the brightness distribu¬ 
tion of the models. Eor example, when selecting rout/^n ~ 25 the 
cut-off occurs at temperatures of ~ 400 K, which corresponds to 
wavelengths around 8 pm. It is therefore important to select rout 
large enough for ah cases to avoid artihcial cut-off in the bright¬ 
ness distribution ( [Honig & Kishimoto[2010[ ). We assume that the 
dust density in the torus declines with distance from the AGN, 
so that at some point a further increase of rout has no impact on 
the SED. We hnd rout/^n ~ 170 is a reasonable choice. This is 
larger than that assumed in most of the previous AGN models. 

The distribution of the dust in AGN is still a matter of 
debate. It was realized early on that the distribution of gas 
and dust is most likely clumpy and in a torus-like structure 
( [Krolik & Begelm^[1988[ ). To account for IR observations 
and the 10 pm silicate feature radiative transfer calcu lations of 
clumpy dust distributions have been carried out, e.g. ([Nenkova 


et al.|200^[Schartmann et al.|2008t [Heymann & Siebenmorgen 

2012| [Stalevski et al.[[2012a[ ). D ust models that assume homo- 

geneous or smooth distrib utions ([Pier & Krohk[1992t[Efstathiou| 
[& Rowan-Robinson|1995|[Efstathiou et al.|2014[ ) have also been 

considered. AGN models assuming clumpy or smooth dust dis¬ 
tributions can produce similar SEDs ( [Eeltre et al.|2012| ). 

Interferometric observations of nearby AGN have so far re¬ 
solved disk-like or extended dust emission components (usually 
coincident with the ionization cones) down to pc scale. Unre¬ 
solved structures such as clumps are below the accessible reso¬ 
lution limits of ~ 15 mas ( [Tristram et al.|2014] ). X-ray monitor¬ 
ing of discrete absorption eclipse events of Seyferts is presented 
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by [Markowitz et alT] ( |2014| ). They find a higher detection prob¬ 
ability of clumps for type 2 than for type 1 , and from the ob¬ 
served column density an optical depth of the clumps between 
Tv = 20 - 16(Qis estimated. The current data of the X-ray study 
provides insight into the cloud distribution only close to the dust 
evaporation zone. For type 2 objects the X-ray absorption can 
also be explained by a homogeneous medium. 

We consider a dust structure in the AGN which is approx¬ 
imated by an isothermal disk that is embedded in a clumpy 
medium. The midplane of the disk is at z = 0, so in the xy-plane 
of the coordinate system the dust density distribution of the disk 
is given by: 


p{x,y,z) =Po — e ( 3 ) 

r 

We choose po so that the optical depth of the disk mid¬ 
plane in the V band, which is measured between rin and rout, 
is Tv,mid = 0,30,100, 300 or 10 00. The scale height is approxi¬ 
mated by ( jPascucci et al.|2004) ): 


8 \rout/ 


(4) 


In addition there are a number Nc\ of clouds that are pas¬ 
sively heated from outside. The possibility that massive stars 
heat clouds from inside is not treated. The density of the clouds 
is assumed to be constant. This is consistent with the small scat¬ 
ter in the column density of the clouds as derived from X-ray 
observations ( [Markowitz et al.|2Q14| ). The clouds are arbitrarily 
distributed in the model space. We assume, for random numbers 
denoted by 3 , that the clouds have an isotropic distribution along 
the azimuth angle 0 = 7 r /2 3 , but they are weighted along the 
polar angle 6 by cos(^) = 3 ^, so that there are more clumps at 
low than at high latitudes, and along the radial direction r an 
equidistant separation with the disc midplane given by: 



A schematic of the AGN model is shown in Fig. The po¬ 
tential accretion disk in the immediate black-hole environment, 
where dust cannot survive, and the jet emission is not treated in 
our model. The dust disk extends from a radius of less than 1 up 
to a few hundred pc distance from the primary heating source. 
We assume dust clouds have a cube size of 0.3 - 3 pc. The hot, 
warm and cold side of these clumps and the shadowing caused 
by them is computed self-consistently in the 3d radiative transfer 
model. Please note that all scales correspond to an AGN lumi¬ 
nosity of 10^^ Lq. As discussed above all radii scale as 


2.3. Technical set-up 

The model grid is set up into (51,51,51) basic cubes with a 
side length of a cube d depending on the choice of the inner 
dust radius. It is given by J = 3.3 rin. The central cube has 
(900,900,900) subcubes, and along each direction the follow¬ 
ing (10,10,10) basic cubes are each subdivided into (15,15,15) 
subcubes. We can achieve a resolution of 4 x 10^"^ cm in the in¬ 
ner part of the torus, whereas by comparison the resolution in 
the hydrodynamical simulations of [Schartmann et al^p014[ ) was 
4 X 10^^ cm. Each cloud is divided into (3, 3, 3) basic cubes and 


Unless otherwise stated the optical depth is specified for the V band. 
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Fig. 2. Visualization of the dust density distribution of an octant of the 
AGN model space. A view from high latitudes (top) and edge-on along 
the x,z axis (bottom) is shown. Clumps deep inside the disk are not 
visible. Note the presence of clumps in the polar region. The axes are in 
arbitrary units. 


each of these is divided into (15,15,15) sub-cubes. This gives 
a total of more than 90,000 cells for each clump, which ensures 
that even for the highest optical depth cases each sub-cube of 
the clump has an optical depth along its axis of tv ^ 1. This 
criterion is similar to the one used by [Honig et~ar ( [2006 ) who 
demonstrated that the cell temperatures converge in that case 
(see their Table l). [Stalevski et al.[ ( [2012a[ ) used clouds with an 
optical depth Ty > 100. These clumps are divided into ( 8 , 8 , 8 ) 
cells when the cloud number is low whereas only a single cell 
is used when the number of clouds is high. The latter simplifica¬ 
tion is not justified whenever the optical depth of a sub-cube of a 
cloud becomes larger than 1. SEDs computed for optically thick 
clouds and ignoring sufficient sub-sampling drastically underes¬ 
timate the IR emission. 


The AGN spectrum (Eq. is divided into 256 frequency 
bins and for each bin 200,000 photon packets are emitted. In 
our treatment this is sufficient to achieve a reasonable signal- 
to-noise of fluxes in the submillimeter. SEDs are computed in 
different viewing angles that are equally spaced in cos 6 between 
0® and 90®. The SEDs are reasonably well sampled by a division 
into nine different viewing angles. 
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2.4. Library setup 

We find that for seeting up a SED library of AGN there are in our 
prescription besides the viewing angle four key parameters: (1) 
the inner radius rin, (2) the volume filling factor of the clouds rj, 
(3) the optical depth of the clouds tv,c 1 , and (4) the optical depth 
of the disk midplane Tv,mid- There are other parameters, such as 
the outer radius, rout = 170 x rtn, up to where dust exists, details 
of the clump structure (Sect. |2.2| ) or their distribution function 
(Eq.|g, which are all kept constant. The variations of the free 
parameters of the SED library are summarized in Table In the 
AGN library there are 400 models each providing SEDs for nine 
viewing directions. 

We varied the density profile within a cloud, e.g. following 
a 1/r profile, and experimented with the cloud structure, for ex¬ 
ample using spherical clumps instead of cubes. We confirm the 
results of |Honig & Kishimoto| ( |201Q| ), who concluded that un¬ 
less very peculiar clump surfaces are used, the exact shape of 
the clumps does not matter too much and a spherical or cubic 
representation probably catches the essence of the dust clouds in 
AGN tori. These authors find somewhat redder MIR colors when 
distributing more dust clouds at larger distances. We also varied 
the cloud distribution along the radial direction and the sizes of 
the clouds. We also experimented with using a large number of 
small clouds in the inner region and a small number of large 
clouds in the outer region. We compared models of the same to¬ 
tal cloud mass and found that all of these changes of the cloud 
properties and distribution functions has some but no significant 
impact on the SED. We therefore do not consider variation of 
these properties in the library. 

Einally, we implemented also clumps with parameters as de¬ 
rived from X-ray eclipse observations ( [Markowitz et al.|2Q14| ). 
Such clumps are detected in the central region of the AGN, so 
in our model in the innermost cube. They are assumed to be as 
small as 3.3 X 10^“^ cm and have ry = 20 - 160. Eor these X-ray 
clumps we computed a single clump temperature. We find that 
clumps detected by X-ray eclipse events have marginal impact 
on the SED as long as their number stays below several million. 
Therefore we neglect such X-ray clumps in the rest of the paper. 


Table 1. Parameters of the AGN library. 


Parameter 

Symbol 

Values 

Viewing angle^^ 

en 

86 , 80, 73, 67, 60, 

52, 43, 33, 19 

Inner radius^^ 

rin (10‘ 'cm) 

3, 5.1, 7.7, 10, 15.5 

Cloud volume 
filling factor^^ 

ri{%) 

1.5, 7.7, 38.5,77.7 

Cloud optical depth^^ 

^V,cl 

0, 4.5, 13.5,45 

Optical depth of 
disk midplane^^ 

^ V,mid 

0, 30, 100, 300, 1000 


Notes. The viewing angle is equally binned in cos(^) and 
measured from the z-axis to the x,y plane. The luminosity 
of the primary source is kept constant at Ln = 10^^L©. The 
cloud volume fraction corresponds to the number of clouds of 
Nc\ = 160, 800, 4000, 8000 within the 3d model space. The 
cloud optical depth is measured along the edge of the cloud that 
has a structure of a cube. The disk midplane is located in the 
x,y-plane at z = 0. We measure Tv,mid along the midplane from 
the source to the outer edge of the disk. 


2.5. Visualization 

As an aid in visualizing the AGN structure we show the dust den¬ 
sity distribution and images at different wavelengths. The partic¬ 
ular parameters of the AGN model we visualize are rin = 0.25 pc, 
ty,c\ = 4.5, Tv,mid = 1000, T] = 1.5 %. A 3d visualization of the 
dust density distribution for an octant of the model space as seen 
from high latitudes and in edge-on view is shown in Eig.[^ The 
density inside the disk is so high that clumps are only visible 
above or within the upper surface layer of the disk. 

Images of the same AGN model are shown in Eig.j^for dif¬ 
ferent viewing directions and wavelengths. The image at 1 pm 
is dominated by scattered light with an exception in the face-on 
view where a small contribution by dust emission is visible in 
the innermost region. In scattered light only the upper and lower 
surfaces of the disk are visible. The optical depth is too high for 
the photons to penetrate deeper into the disk. This explains why 
there are two distinct gray areas noticeable. Clumps are visible 
as the few bright dots in the middle and bottom panel (left). The 
image at 10pm displays emission by warm dust and scattering is 
no longer important. The optical depth becomes smaller and one 
probes deeper inside the disk (compare top-left and top-middle 
panel of Eig.[^. The disk itself remains optically thick at lOyum 
so that two distinct surfaces are visible as in the scattering light 
images. As can be seen in the top-middle panel, the emission 
from the diffuse dust in the ionization cones can dominate the 
total emission from an edge-on view. This can explain why some 
Seyferts show extended emission coincident with the ionization 
cones at 10pm despite the fact that there is an optically thick 
torus which lies in a plane perpendicular to the ionization cones. 
The image at 10pm shows also the shadowing of the clumps. 
They can be noticed as notched contours or jet-like structures 
(top-middle). Behind a clump there is less direct heating from 
the central source, the dust becomes colder and the lOyum flux 
is reduced. The image at 100yum reveals the cold dust emission 
and at this wavelength the disk is optically thin and the radiation 
appears isotropic. 

The size of the torus becomes wavelength dependent if 
one cuts the image contras t at higher flux level s, e.g. take 
log(E/Eniax) > -4 in Eig. [3] ( [Stalevski et al.|2012b| ). At shorter 
wavelengths radiation from the inner region dominates and at 
longer wavelengths, the emission arises from cold dust deeper 
inside the disk and from further out. 


2.6. AGN Dust Properties 


Dust in dense environments such as in the circumnuclear re¬ 
gions of AGN is certainly different from the dust in the diffuse 
ISM. The variation is caused by the huge difference in the den¬ 
sity in both media and for AGN some observational evidence of 
a chan ge in the dust properties is presented by [Maiolino et al. 
( |2QQ1| ). The basic process for this modification of the dust is 
probably grain coagulation which leads to fluffy particles and 
grain growth ( [Mathis & Whiffen|1989| ). 

The 2175 A extinction bump, which is attributed to small 
graphites and PAHs, is not detected in AGN ( [Maiolino et al. 
[200 1[ ). The lack of PAH emission in the proximity of the AGN 
was used as an indication that particles smaller than lOnm are 
destroyed ([Voit[[1992t [Mason et al.[[2007[). However, [Alonso- 
[Herrero et al.[ ( [2014] rand [Ramos Almeida et al.[ ( |2014b[ ) detected 
in 6 Seyfert galaxies the 11.3 pm PAH feature in the nuclear re¬ 
gion. They show that PAH molecules can survive in the nuclear 
environments of low luminosity AGN as close as 10 pc from the 
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Fig. 3. Images of the AGN torus computed by ray-tracing at viewing directions corresponding to 15° (face-on, bottom), 45° (middle) and 75° 
(edge-on, top) and from left to right at wavelengths of 1 (grey), 10 (orange), and 100yum (blue). The images are logarithmically scaled and at each 
wavelength the flux is normalised to its peak value. 


central source. They interpret this as evidence that the PAHs are 
not destroyed in the torus by the hard radiation of the AGN. 


Inner radii derived from K-band reverberation mapping of 
type 1 AGN are a facto r of 3 smaller than wha t is expected 
from standard ISM dust. |Kishimoto et al.| ( |2QQ7| ) suggest that 
this is due to the fact that the grains are larger than standard 
ISM dust. However, |Honig & Kishimot^ ( |2010| ) conclude that 
the grain size distribution in AGN is the same as that usually 
adopted for the ISM. Previous models of AGN tori assumed dif¬ 
fuse ISM dust, graphite dominated dust, and dust with different 
optical constants than usually assumed for ISM dust. Fritz et al. 
(2006) shows that variations of the silicate feature strength may 
be explained by changing the optical constants of the dust or its 
distribution in the torus. 


In this study we consider fluffy mixtures of silicate and amor¬ 
phous carbon grains ( [Krugel & Siebenmorgen||1994|). Op tical 
constants of the materials are from [Zubko et al.| ( |1996| ) and 
|Draine| ( [2003| ), respectively. The bulk density of the dust mate¬ 
rials is 2.5 g/cm^, the carbon-to-silicate abundance ratio is 6.5, 
and 50% of the volume fraction of fluffy particles is vacuum. 
We consider large particles with radii a between 16 and 260 nm. 
They follow a grain size distribution with number density oc 
( [Mathis et al.|1977| ). The influence of the dust cross section on 
the AGN emission is exemplified in Fig. For comparison we 
apply dust cross sections of pure (unmixed) particles that we call 
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ISM dust ( jSiebenmorgen et J^|2014| ). Besides the grain struc¬ 
ture, ISM grains are otherwise identical to fluffy grains. They 
have the same optical constants, abundances, and sizes. Fluffy 
grains are more efficient absorbers and have lower scattering 
cross sections than ISM dust for the whole spectrum. The extinc¬ 
tion cross section of fluffy grains in the V band is = 33.3 (g- 
dust/cm^). This is a factor of 2 larger than for ISM dust. At 1 mm 
the cross section of fluffy grains is 12 times that of ISM dust. 
For the purpose of a comparison we take an AGN library model, 
computed using fluffy grains, inner radius rin = 0.25 pc, cloud 
filling factor p = 1.5%, cloud optical depth tv,c 1 = 4.5, and 
optical depth of the disk midplane of Tv,mid = 1000. A model 
with the same dust density distribution is run assuming ISM 
dust. For the latter model the disk midplane optical depth re¬ 
duces to Tv,mid = 750 and the optical depth of the clouds to 
^v,c\ = 3.375. Fluffy grains emit more strongly in the far IR 
and submillimeter than ISM dust and the SED peak wavelength 
shifts towards longer wavelengths. The cross-over wavelength 
where the AGN emission in face-on and edge-on views is iden¬ 
tical, in other words the wavelength where the torus emission 
becomes isotropic, occurs for fluffy grains in the submillimeter 
(~ 200yum) and for ISM dust in the far IR (~ 100yum). 

In face-on view the emission in the mid- and near IR is sim¬ 
ilar whereas in edge-on view the small difference in the extinc¬ 
tion of both dust models leads to a significant change in the con- 
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Fig. 4. Extinction (full line) and scattering (dashed) cross section K 
as a function of wavelength {top). We show fluffy aggregates (bl ack) 
of silicates and amorphous carbon ( [Krugel & Siebenmorgen|19^ and 
identical materials as unmi xed bulk particles (red) of the diffuse ISM 
( [Siebenmorgen et al.|2014| ). Eor both dust models and with otherwise 
identical parameters the emission spectrum of a clumpy AGN torus is 
given for face-on and edge-on view {bottom). 


tinuum flux. In the silicate band at lOyum the fluffy grains have a 
shallower and broader absorption cross section when compared 
to ISM dust. This is also reflected in the SED that is shown in 
Fig. 1^ In edge-on view the silicate absorption is noticed using 
ISM dust and not when fluffy grains are considered. In face-on 
view the silicate feature has a peak emission at 10.4 pm for ISM 
dust and at 11.5 pm for fluffy grains. In the dense dust environ¬ 
ment of the AGN the particles have a higher sticking probability 
than in the diffuse ISM and form rather fluffy dust agglomer¬ 
ates than dust in the diffuse ISM. The 2175 A extinction bump is 
not observed in AGN spectra and therefore we use in both dust 
models amorphous carbon that does not produce such a feature. 
In the scattering cross section ISM dust shows a feature at 0.2 pm 
( Hofmeister et al.|2009| ) that is not present for fluffy grains. Such 
a scattering peak is not detected for AGN. 


2. 7. AGN Extinction 

The fact that dust near AGN is different from the diffuse ISM 
may be best witnessed from changes in the extinction curve. 
However, we point out some observational difficulties. For un¬ 
resolved observations, photon scattering within the beam may 
increase the detected flux altering the wavelength dependence 
of the extinction, and for clumpy media there are light paths to¬ 
wards the source at different optical depths. These effects are 
known to change the appearance of extinction curves and to ac¬ 


count for this the effective o ptical depth is introduced ( jNatta & 
|Panagia| 19841 |Krugel|2009| ) : 


Teff ■ 


In 


^obs 

Fnd 


( 6 ) 


where Fobs is the observed flux and Fnd is the flux that would 
be observed in the absence of dust. Both quantities are easily 
derived by counting the number of photons within the aper¬ 
ture of an AGN model and comparing them to an identical MC 
run but w ithout dust. We illustrate this for both dust models of 
Sect. 12.61 and take as the AGN model: r^n = 0.25 pc, tv,c 1 = 4.5, 
TY,mid = 1000, T] = 1.5 %, and in addition a disk-only model 
without clumps {p = 0). For these four models it was neces¬ 
sary to improve the photon statistics in the optical. We do this by 
launching 40 times more packets than for the SEDs computed 
for the library. 


The effective extinction curves are displayed in Fig. for 
three directions: face-on {6 ~ 33°), an intermediate case with 
^ ~ 50°, and edge-on {6 > 66°). For reference we also show the 
extinction curve as derived from the dust cross section F/Fy of 
the fluffy grains. In pencil-beam observations the effective ex¬ 
tinction and the dust extinction are identical. It is striking that 
none of the effective extinction curves resembles the one of the 
dust model. There is no direct one-to-one link between the ob¬ 
served extinction curve and the wavelength dependence of the 
dust cross sections (Scicluna & Siebenmorgen 2015, subm.). 


In edge-on view Tv,eff > 0 from the UV down to the near IR. 
In the UV the extinction curve flattens, stays below the extinction 
of the dust model and becomes eve n gr ay in accordance w ith the 
findings by jNatta & Panagia| ( |1984| ) and jWolf et al.| ( |199^ . From 
similar observational results it was claimed that the dust compo¬ 
sition in the circumnuclear region of AGNs could be dominated 
by large grains ( jMaiolino et al.pOOT] ). As demonstrated in Fig.[^ 
such a statement does not necessary hold. One cannot directly 
conclude that there is an increase of grain size when observing a 
flatter than usual extinction curve and this is especially the case 
when the object remains highly unresolved. 


The effective extinction may also become negative. This oc¬ 
curs whenever there is more light detected with dust than without 
dust ( |Krugel|2009| ). Negative extinction can be measured when 
extra light is scattered into the beam and overshines the primary 
source. For intermediate or face-on views presented in Fig. 
this occurs at wavelengths >0.8 yum. In this case there is an op¬ 
tically thin view to the source and in addition light is scattered 
from the disk into the beam. In the face-on view of the disk-only 
model the enhanced scattering contribution is important already 
in the V band so that Tv,eff < 0. This explains the behaviour 
of the dashed red line in Fig. and remains stable as long as 
0 < 44° (close to face-on). For edge-on views variations in 6 
are unimportant as long as the disk provides substantial extinc¬ 
tion. For clumpy models and viewing angles that do not hit the 
disk {0 < 60°) small variations in 6 have a large impact on Tv,eff, 
because by chance a clump may obscure the source or not. 


The effects described above become more pronounced for 
ISM dust (dotted) than for the fluffy grains because the former 
have stronger scattering efficiencies (Fig.|^ top). They also pro¬ 
duce significant steeper extinction curves than is expected from 
the dust model, in agreement with the qualitative estimates de¬ 


rived by [Krugel] ( |2009| ). 
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Fig. 5. Impact of AGN parameters on the SED in face-on (top) and edge-on (bottom) view. We vary the inner radius rin for models a) without and 
b) with a disk, c) the optical depth of the disk midplane Tv,mid, d) the volume fraction of the clouds t], and e) the optical depth of the clouds tv,ci- 
Unless otherwise specified rin = 0.25 pc, Tv,mid = 100, rj - 38.5%, and tv,ci = 4.5. 



Fig. 6. Wavelength dependence of the effective extinction curves for 
three viewing directions: edge-on (black), intermediate (blue) and face- 
on (red) of AGN models as described in the text. The effective extinction 
curves of the AGN model computed with fluffy grains are shown as full 
line and for disk-only models as dashed lines, that of ISM dust as dotted 
lines. The extinction curve (thick gray) derived from the cross section 
of the fluffy grains is given for reference. All curves are normalized to 
the V band. 


3. Impact of parameter variation on the SED 

We illustrate in Fig. [^changes in the SED when one parameter 
of the AGN library is varied while all the others stay fixed. For a 
particular model the SED in face-on view is shown in the top and 
in edge-on view in the bottom panel of Fig.In both views, and 
respective panels, the SEDs are identical when the AGN emits 
isotropically. This happens at wavelengths long-ward of the far 
IR peak flux. The AGN torus emits an-isotropically when the 
SED strongly depends on the viewing direction. This occurs at 
wavelengths short-ward of the far IR peak. Panels a) and b) of 
Fig. [^illustrate how an increase of the inner radius shifts the far 
IR peak to longer wavelengths. When the inner radius becomes 
larger the dust gets colder, hence shows a maximum emission at 


longer wavelengths so that in the far IR and submillimeter the 
flux increases while in the near IR the flux decreases. 

Clumpy models without and with a homogeneous disk are 
shown in Fig.[^a,b, respectively. In face-on view models with a 
disk show in the near IR a huge flux increase when compared to 
models without a disk. The addition of disk material increases 
the total dust mass and this leads to an increase of the sub¬ 
millimeter flux at 200yum by a factor 2-3. The impact of the 
midplane optical depth of the disk Tv,mid is shown in Fig. [^c. 
For edge-on views the disk provides additional dust extinction 
changing the silicate band from emission to absorption, while 
for face-on views the influence of the disk on the strength of that 
feature is small. When disk material is included there is addi¬ 
tional hot dust in the system, and in face-on views the near IR 
(~ 2 yum) flux increases by orders of magnitude (Fig. [^c, top). 
The near IR flux may increase even further when a puffed-up in¬ 
ner rim, as known to exist in proto-planetary disks, is considered 
( [Dullemond et al.|2002[|Siebenmorgen & Heymann|201^ . The 

disk reduces the near IR flux for edge-on views provided that the 
extinction is high enough (Tv,mid > 30, Fig.[^c, bottom). The op¬ 
tical depth of the clouds tv,c 1 also has a huge impact on the near 
IR and submillimeter fluxes (Fig.j^e). In face-on view a change 
of the silicate band from emission to absorption is noticed when 
the optical depth of the clouds tv,c 1 is increased (Fig. [^e). For 
edge-on views an increase of p will increase the dust extinction 
and hence produce a stronger 10 yum absorption feature (Fig. |5]d, 
bottom). The impact of par ame ter variation on the silicate band 
is further discussed in Sect. 13.11 

3.1. Silicate Feature 


Mid IR spectra of AGN show silicate features in absorption 


(|Roche et al. 19911 |Spoon et al. 2002[ |Siebenmorgen et al. 

2004bl |Levenson et al. 2007|) oi 

' in emission (|Siebenmorgen 

et al.|2005[ |Hao et al.||2005| 20071 Sturm et al.||2005l |Sirocky 

et al.|2008 ).|Hatziminaoglou et al. 

(|20151) present a census of the 


silicate features in mid IR spectra of 800 AGN observed with 
Spitzer. Where the feature is in emission, in about 65% of the 
cases the peak wavelength is > 10.2yum whereas the shift of 
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Fig. 7. Histograms of the silicate emission (red) and absorption (grey) band for the AGN library labeled “Disk-rClouds” (left) and disk-only 
models (right). Distribution functions are shown from top to bottom for the the central wavelength of the feature, the viewing angle, the inner 
radius of the AGN dust torus, and the optical depth along the midplane of the homogeneous disk . 


the 9.7yum absorption features is much smaller. The 10/im sili 
cate emi ssion and absorption features are also observed from the 
groun d ([Roche et J^|2006 Honig et al.||2010l lAlonso-Herrero 


et al. [201 2[ [Sales et al.||2011[ [Esquej et al.||2()14[ [Gonzalez- 

Martin et al.|2013t|Ramos Almeida et al.|2014bt|Ruschel-Dutra 


et al.||2014| ), and it is noticed that the depth of the lOyum sili¬ 

cate absorption band is enhanced compared to that measured by 
Spitzer. 

In the MC code we applied a fine wavelength sampling that 
allows us to discuss the behaviour of the lOyum band. The dust 
cross section of the lOyum silicate band is larger than that of the 
18yum feature (Fig. |^. The feature strength can be quantified 
similarly to the definition of Eq. ([^ by its effective optical depth 


Tsi = - In 


-^peak 

-^cont 


(7) 


where Fpeak is the maximum or minimum of the observed in- 
band flux at center wavelength and Fcont is an estimate of the 
underlying continuum at that wavelength. We estimate the con¬ 
tinuum for each model spectrum by a straight line with anchors 


set outside the red and blue wing of the band at wavelengths 
between 5 - 7yum and 12 - 14yum. In this definition, emission 
features appear at rsi < 0 and absorption features at rsi > 0 . 
Clumpy media offer some direct views to the inner hot regions 
that fill the absorption and reduce the depth of the band. The fea¬ 
ture strength depends on the distribution of the dust within the 
AGN torus ( [Levenson et al.|2007] ), the dust mineralogy ( [Sirocky 
|et al.||20()8] ), the grain structure (Fig. |^, and is rather insensi¬ 
tive to the source spectrum. Note that, similarly to the discussion 
of the extinction curve (Sect. |2.7| ), the effective optical depth de¬ 
rived from the 10 fim silicate absorption feature does not provide 
a measure, and at best is only a lower limit, of the optical depth 
of the dust. 

Histograms of the strength of the lOyum silicate band are 
shown in Fig. for the AGN library as a whole and disk-only 
models. The latter are selected from the library of models with¬ 
out clumps (77 = 0 ). 

For Seyfert Is the observed feature strength ranges for emis¬ 
sion bands between -0.15 < rsi < -0.01 and for absorp¬ 
tion bands between 0.06 < rsi < 1.5, for quasars between 
-0.3 < Tsi < -0.01 and rsi < 0.62, respectively ( [Thompson 
et al.||2009 ). Schartmann et al. ( 2008| ) observe a strength of the 
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Fig. 8. Distribution of the strength of the 10 pm silicate band (Eq. for 
the AGN library (labeled “Disk-rClouds”, gray) and for the disk-only 
models (red). Emission features have Tsi < 0 and absorption features 
Tsi > 0. The low frequency tail of the distribution function is truncated 
atTsi = 1.4 



absorption band for type 2 AGN between 0.2 < rsi < 3.5. The 
AGN library accounts well for the observed scatter (Fig. [^. In 
the setup of disk-only models absorption features stronger than 
Tsi >0.8 are not present. [Levenson et al.| ( |2007] ) report a mean 

feature strength for Seyfert Is around tsi-0.2 and for Seyfert 

2s around fsi ~ 0.2. Average values derived from the distribu¬ 
tions functions presented in Fig.j^are about the same. We cannot 
identify a striking need to postulate clumpy AGN models when 
explaining observed feature strengths of the 10 pm silicate band. 
Arbitrary clump distributions ease accounting for silicate fea¬ 
tures observed at any strength and viewing angle of the AGN. 

Disk configurations where the 10 pm band is in absorption 
or emission can be seen in the distribution functions presented 
in Fig. [7] The distribution of the central wavelength position 
Ac of the silicate band is shown in the top panels of Fig. □ in 
the silicate band the cross section ^si peaks at 9.5 pm and the 
exact position depends on the choice of optical constants. Ab¬ 
sorption features detected in the SED match this position with 
some scatter that is explained by radiative transfer effects. No¬ 
ticeably for most models Ac of the emission feature is shifted 
to 10.5 - 11.9yum. The shift of the extrema to such long wave¬ 
lengths is naturally explained assuming that the emission fea¬ 
ture is optically thin, so that the observed flux is estimated by 
F ^ K X B(T^). The folding of K with the steeply rising Planck 
function B(T(^) of the dust at temperature results in the ob¬ 
served wavelength shift (|Siebenmorgen et al.|2005[|Nikutta et al.| 
|2009l ). 


The silicate feature may appear in emission whenever there 
is an optically thin view to the inner, hot regions of the disk. This 
may happen for any Tv,mid of the disk midplane, inner radii and 
viewing directions (Fig. right). Most of the emission bands 
in disk-only models are seen close to face-on views. There are 
only a few disk-only models that have sufficient dust extinction 
so that the silicate band is in absorption. These are models with 
Tin > pc/4. Tv,mid ~ 100, and the disk view needs to be ^ > 60®. 
When clumps are considered (Fig. [7] left) the silicate band may 
appear either in absorption or emission for an almost arbitrary 
choice of disk configurations and viewing angles. It is striking 
that there are clumpy models displaying the emission band in 
edge-on view. There is a rather constant decrease of the number 
of ’Disk-i-Clouds’ models that show the emission band when the 
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inner radius or the optical depth is increased or vice versa for 
the absorption band. All this is consistent and explained by the 
corresponding increase or decrease of the dust extinction and 
temperature. 


3.2. Apparent versus intrinsic AGN luminosity 

Because of the assumed spectral shape of the AGN (Eq.|^, pho¬ 
tons are predominantly emitted from the central source at far UV 
to near IR wavelengths. Dust near the AGN scatters and absorbs 
these photons and re-emits this radiation in the infrared. By inte¬ 
grating over the observed fluxes one derives the ’apparent’ AGN 
luminosity Lobs- One is however interested in deriving from ob¬ 
servations the intrinsic luminosity of the AGN source Lagn- The 
AGN library provides a utility that allows the conversion of the 
apparent (observed) luminosity to the intrinsic luminosity of the 
AGN. 

In the models the apparent luminosity is computed from pho¬ 
tons emitted into a given viewing bin. We compute the observed 
luminosity for models where dust is considered or neglected. For 
the latter, the no-dust models, the luminosity is constant because 
the central source emits isotropically. Therefore in one of the 
nine viewing bins used, the luminosity Lii/9 = Lnd(^), 

where in the library D = 50Mpc, and F^^{0) is the flux inte¬ 
grated over frequency in a particular 6 direction for the “no dust” 
model. For models with dust the luminosity is computed from 
the flux as given in the library, and Lobs(^) = 47rD^ F(6). We 
define the anisotropy conversion factor 



( 8 ) 


The conversion factor becomes larger than 1 whenever there 
are more photons detected in a viewing bin than originally emit¬ 
ted in that direction from the central source. The conversion fac¬ 
tor 6 allows one to derive the intrinsic AGN luminosity Lagn- 
The ratio of the observed-to-intrinsic luminosity is given by 


Tagn - 9 


Tobs(^) 


(9) 


The distribution functions of Lobs/L agn are shown in Fig.[^ 
for different AGN models. Typically the apparent luminosity is 
about two times larger for type Is (face-on) than for type IIs 
(edge-on). The fact that for the same intrinsic AGN luminos¬ 
ity, models which correspond to type 1 sources emit more IR 
radiation than type 2s is clearly visible in the SEDs (Fig. [^. 
Such an increase of the IR radiation is observed in a sample of 
SCR sources ( [Siebenmorgen et ar]|2004a| ). For face-on views 
there are often more photons radiated into the beam than orig¬ 
inally emitted from the source, hence Lobs/L agn > 1- Models 
with clumps show for edge-on views a broader distribution in 
Lobs/L agn than models without clumps. The impact of the cloud 
filling factor j] is shown for models without a disk in Fig. [^(bot¬ 
tom). For viewing directions with a single cloud along the line 
of sight typically 40% of the AGN luminosity is observed. This 
may happen for AGN without a disk and with a small number 
of clouds (rj = 1.5 %). When there are more clouds Lobs/LAGN 
increases and ranges between 0.7 - 1.3. 

We note that the ’apparent’ AGN luminosity, which is usu¬ 
ally what is quoted in the literature, must simply be divided by e 
in order to estimate the intrinsic luminosity of the AGN Eq. (^. 
The conversion factor e is the inverse of the anisotropy factor A 
defined by |Ef stathiou| ( [2006] ) . 
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Fig. 9. Histograms of the apparent-to-intrinsic AGN luminosity (Eq.[^. Results are given for the complete AGN library (top), and for disk-only 
models (middle) in face-on (6 ^ 43°) and edge-on (6 ^ 73°) views. The luminosity ratio is shown for the volume filling factor of the clouds rj 
(bottom) for models without the homogeneous disk. 


4. Testing the SED library 

In this section we fit the observed SEDs of a number of repre¬ 
sentative objects from the near IR to the submillimeter including 
Spitzer/IRS spectroscopy, Herschel photometry of distant high 
luminosity AGN and high spatial resolution NIR and MIR ob¬ 
servations of nearby, low luminosity Seyfert nuclei. Further we 
test that sources of the library would be classified as AGN when 
applyi ng the empirical mid IR selection criteria by [Stern et aL] 
( |2Q05| ). First, we test how much of the observed SED can be ex¬ 
plained by pure AGN activity ( [Fritz et ar]|2006| ). As a second 
step, we combine the predicted SED from the AGN library with 
an additional starburst (SB) component. We call such combina¬ 
tion ’AGN-fSB’ models. To this end we select elements from 
the SE D library of starburst models of [Siebenmorgen & Krugel[ 
( 2007[ ). We consider the starbursts radius Rsb, the luminosity 
Lsb, the visual extinction Ay, and the ratio of the luminosity that 
is due to OB stars t/ob as free parameter and keep the hot spot 
density constant t/hs = 10"^cm“^. In addition, to better match 
data below 2 /jm we add an extra component. For quasars we 
use beamed synchrotron emission, which is approximated by a 
power-law, and for the galaxies stellar light from the host, which 
is approximated by a blackbody. The extra components add a mi¬ 
nor contribution to the photometry between the thermal IR and 


the submihimeter. Nevertheless we subtract the flux of the extra 
components from the data when fitting the dust emission. For the 
latter fit we use a Bayesian for malism that utilizes t he Markov 
Chain Monte Carlo algorithm by [Johnson et al.[ ( [2013[ ). This gen¬ 
eral purpose SED fitting tool employs the Metropolis-Hastings 
algorithm, e.g. [Metropolis et aL] ( [1953| ). The code is able to take 
any set of SED templates. In our case we use either the pure 
AGN library or simultaneously the combination of the AGN-hSB 
libraries. The tool returns the requested best-fitting parameter es¬ 
timates, and in addition construct the confidence levels from the 
posterior parameter distribution. We verified the accuracy and 
robustness of that SED fitter against resul ts obtained with tradi- 
tional least-squares methods. We confirm [Tohnson et'ar]p013| 
that the programme recovers best-fitting values similar to those 
from traditional methods. 

In the fit procedure the Spitzer spectroscopy data are binned 
to a spectral resolution that adequately reproduces the shape of 
the SED including the silicate feature. We display the SED in 
Figs. [^-[^ and the best fit parameters including their uncer¬ 
tainties are summarised for the pure AGN models in Tableand 
for the combination of AGN and starburst activity in Table 
Apart from the inclination, the fit parameters may cover almost 
the entire parameter range because of the arbitrary cloud dis¬ 
tribution. One finds adequate fits (solutions) of the SED using 
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Fig. 10. Re st frame SEP of the quasars 3C249.1 (top) and 3C351 (bottom). Data are from the NASA/IPAC Extragalactic Database (NED), 
Spitzer/IRS ( [Leipski et al.|2010] ) and Herschel (this work. Table |^. A fit using a pure AGN is shown on the left and for a combination (black line) 
of AGN (red line) and starburst (blue dashed) activity on the right. A synchrotron component (yellow dots) is added to better match the near IR. 
Model parameters are given in Table [^. 


either a combination of a large number of clouds with small op¬ 
tical depth, or models that have a small number of clouds with 
high optical depth. Such a change of the cloud distribution re¬ 
quires that the inner radius, the optical depth of the disk or other 
of the free parameters are modified as well. 


ble[^ The total flux is measured and uncertainties are computed 
by adding on the statistical error a calibration uncertainty of 7% 
for PACS and 5.5% for SPIRE bands as well as a confusion noise 
of 5.8, 6.3, and 6.8 mJy in the bands at 250, 350, and 500yum, 
respectively. 


4.1. Observations 

We derive aperture photometry on the flnal images extracted 
from the Herschel archive and apply correction factors for the 
encircled energy fraction and colors as given by [Balog et~ar| 
( 2014| ) and the User Manuals. The background is subtracted and 
estimated from the mean flux that is measured in a 2" wide an¬ 
nulus which is placed 1" further out than the source aperture. 

PACS and SPIRE photometry of 3C249.1 and 3C351 (PI: 
Luis Ho) is available in 6 Herschel Alters and summarised in Ta- 


PACS photometry of the Seyfert nuclei NGC 1365 (PI: M. 
Sanchez), NGC 3081 (PI: M. Sanchez), and NGC 4151 (PI: A. 
Alonsoh) is measured by adjusting a radius of the source aper¬ 
ture that matches the diffraction limit of the Herschel telescope. 
We apply so urce radii of 6", 1", and 12" at 70, 100, and 160yum, 
respectively ( [Balog et"ST]|2014| ). We centre this aperture on the 
location of the AGN. These galaxies are resolved by Herschel 
but their AGN tori are not. Eor the other objects Herschel pho¬ 
tometry is obtained from [Efstathiou et al.[ p013| ) and [Podiga- 
[choski et'ar] ( |2015[ ). 
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Fig. 11. Rest frame SED o f the radio galaxies 3C324 (top) and 3C356 (bottom). Data are from NED, Spitzer/IRS ( [Leipski et al.|2010| ) and Herschel 
( [Podigachoski et al.|2015| ). A fit using a pure AGN is shown on the left and for a combination (black line) of AGN (red line) and starburst (blue 
dashed) activity on the right. Stellar light from the host (yellow dots) is added to better match the near IR. Model parameters are given in Table 


4.2. Fits to the quasars 3C351 and 3C249.1 

The near IR photometry of 3C351 is fit assuming a contribution 
of a beamed synchrotron component that we approximate as a 

power-law (Fy oc y^) with a slope a -0.2. We show that pure 

AGN models with parameters are summarized in Tablej^account 
for the complete SED without the need of postulating additional 
starburst activity. 

The best fit of the AGN+SB models shows that the SED can 
be explained mainly by a torus and a small (8 %, Tablecontri¬ 
bution from a starburst. The predicted viewing direction is close 
to face on (6 ~ 33°). The models reproduce the silicate emission 
feature of this quasar and the near IR bump without the need 
for an additional hot dust component as suggested by |Mor et al.| 
( 2009| ). This is a direct consequence of additional material in the 
innermost region of the torus that is provided by the homoge¬ 
neous disk and by clouds in the ionisation cones of the AGN. 
We note that our conclusion that the starburst makes a negligible 


contribution to the bolometric luminosity is a consequence of the 
large ratio of Tout/^in ~ 170, and the assumed fluffy grain model. 
The far IR emissivity of fluffy grains is much larger than that of 
ISM dust (Sect.[^ Eig.|4]). 

Eor 3C249.1 we And that pure AGN models provide a rea¬ 
sonable fit to the SED and its silicate emission feature without 
the need for SB activity. There are of course AGN-fSB models, 
similar to 3C351, that fit the data. Again, we assume synchrotron 
emission that makes a significant contribution to the near IR. We 
And that in this quasar the contribution to the total luminosity 
from a starburst is 3 % (Table [^. 

4.3. Fits to the radio galaxies 3C324 and 3C356 

The first thing that is striking about the two radio galaxies we 
studied here (3C324 and 3C356) is that their SEDs are distinctly 
different from those of the quasars. They show a much weaker 
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Fig. 12. As Fig.llljfor the hyperluminous infrared galaxy IRAS 10214+4724 with data from NED, Spitzer IRS ( [Hernan-Caballero & Hatzimi-| 
|naoglou|2011| ), and MIPS, IRAC and Herschel by [Efstathiou et al.] ( |2013| ). A fit using pure AGN models are shown with (dashed) and without 
(dotted) clumps. Model parameters are given in Table 


Table 2. Herschel photometry with Icr fiux uncertainty. 


Object 

Eilter 

Oum) 

Flux 

(Jy) 

NGC 1365 

70 

128 + 14 


100 

148 + 16 


160 

179 + 20 

NGC 3081 

70 

2.8+ 0.3 


100 

3.1 +0.4 


160 

3.5+ 0.4 

NGC 4151 

70 

5.2+ 0.6 


100 

3.4 + 0.4 


160 

4.9 + 0.6 


Oum) 

(mJy) 

3C249.1 

70 

87.9 + 6.2 


100 

70.9 + 5.7 


160 

40.6 + 3.3 


250 

30.7 + 8.0 


350 

41.7 + 8.7 


500 

50.7 + 9.4 

3C351 

70 

189.6+19.1 


100 

164.4 + 16.2 


160 

78.5 + 9.3 


250 

48.9 + 9.1 


350 

19.6 + 7.5 


500 

16.5 + 9.2 


near- to mid IR continuum relative to the far IR, a completely 
different slope in the 1 - 10 pm range and evidence of shallow 


absorption features at lOyum. All of these differences are under¬ 
stood in terms of the optically thick torus model and are strongly 
supportive of the idea of unification. For the radio galaxies we 
add to the pure AGN torus or the AGN+SB models contribu¬ 
tions from black bodies at temperature of 4000 - 5000 K, which 
model the emission from old stars in the host galaxies in these 
objects. The SED of both galaxies are fit by pure AGN models 
with parameters as given in Table In the combined models the 
starburst component is 8 % for 3C324 and 19 % for 3C356 of the 
intrinsic AGN luminosity (Table |^. 


4.4. Fit to the hyperluminous infrared galaxy IRAS 
F10214+4724 


We discuss a fit of the SED of the hyperluminous infrared galaxy 
IRAS FI0214 +4724 which has an apparent luminosity that ex¬ 
ceeds 10( Rowan-Robinson et al. 1 1991 ). It is generally ac¬ 
cepted that the enormous luminosity of this galaxy at redshift 
z=2.224 is due to gravitational lensing by a foreground galaxy 
with magnification betwee n 50 and 1 00 ([Broadhurst & Leharj 
|1995t [Serjeant et al.|1995| ) or 15 - 20 ( [Deane et al.|2013|). Fol 


lowing the detection of a silicate emission feat ure ([Teplitz et aT 
2006[) in this narrow -lined object, [Efstathiou] ( [2006[ ) and [Efs¬ 


tathiou et al.[ ( [2013| ) proposed that the emission feature is due 


to emission from narrow-line region clouds at a temperature of 
210K. 


Here we find that the SED of IRAS E10214+4724 can also 
be fit by a pure AGN model without clouds when a small in¬ 
ner radius or in a ’Disk+Clouds’ configuration when larger in¬ 
ner radii are assumed (Table [^. In the AGN+SB models we find 
the ratio of starburst to AGN luminosity is 20 % (Table |^. The 
model presented here for IRAS El0214+4724 is similar to that 
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Fig. 13. SED fits using pure AGN models (magenta) of Seyfert nu¬ 
clei: NGC 1365, NGC 4151, NGC 3081, and NGC 5643. High spatial 


resolution photometry is from Ramos Almeida et al.| ( |2014b|). Ground 
based MIR spectroscopy of NGC 1365 is performed by|Burtscher et al.| 
( |20T3l ), NGC 3081 and NGC 4151 by [Ramos Almeida et al.| ( |2014b| ), 
and NGC 5643 by |Honig et al.| ( |2010| ), respectively. FIR photometry of 
the nuclear region provided by NED or Herschel (Table are plotted 
as upper limits of the AGN torus emission. Model parameters are given 
in Table |3] 


4.5. F/fs to Seyfert nuclei 

|Garcia-Burillo et ar] ( |2014| ) presented ALMA continuum obser¬ 
vations of NGC 1068 and find that rout/^n ~ 100. This is much 
larger than the value inferred by |Alonso-Herrero et~ar] (|2011| ) 
(^out/^in ~ 6 ) who used MIR data. [Ichikawa et al.| ( |2015| ) also 
on the basis of NIR and MIR data, find that this ratio is less than 
20 for a sample Seyferts. It is expected and confirmed by our 
models that the ratio derived from submillimeter data, which is 
probing dust at ~ 30 K or below, will be larger than the ratio 
derived from MIR data which is probing dust at 300 K. 

[Lira et al.] ( |2013| ) find that in 50% of the Seyferts that are n ot 
well fit by the clumpy torus model of jNenkova et al.| ( [2010| ), a 
NIR excess is observed. In our models such a NIR emission is 
provided by the homogeneous disk component. 

Ramos Almeida et al.j ( |2014b| ) present high spatial resolu 


tion 1 - 18yum data of a sample of nearby undisturbed Seyfert 
galaxies with low to moderate amounts of foreground extinc¬ 
tion (Ay < 5 mag). We model four objects from their sample (the 
Seyfert Is NGC 1365, NGC 4151 and the Seyfert 2s NGC 3081 
and NGC 5643) to test the ability of our SED library to fit 
available nuclear NIR and MIR photometry and spectroscopy. 
We take high spatial resolution MIR spectroscopic observations 
of NGC 1365 from [Burtscher et al.j (|20T3l), NGC 3081 and 
NGC 4151 fromjRainos Almeida et al.| ( [M^ , and NGC 5643 


from jHonig et alT (2010). Our results are shown in Fig. p3] and 
the fit parameters are summarized in Table The torus models 
clearly under-predict the far infrared data (plotted as upper limits 
in the figure) as the emission at these wavelengths is dominated 
by the host galaxy. 


4.6. Two color I RAC diagram 

Mid IR photometry provides a simple technique for identifying 
active galaxies. [Stem et aL](|2005|) repo rt on Spitzer mid IR col¬ 
ors, derived from IRAC ( |l4zio ct al.||2004| ), of nearly 10,000 
spectroscopically identified AGN and galaxies. They find that 
simple mid IR color criteria provide remarkably robust sep¬ 
aration of AGN from normal galaxies and stars. We indicate 
the empirical relation in Fig. In that diagram AGN popu¬ 
late the area in between the dashed lines at [3.6] - [4.5] > 0.3, 
whereas galaxies and stars are below at [3.6] - [4.5] < 0.3 and 
-0.3 < [5.8] - [8.0] < 3.2. 

We compute the IRAC colours for each spectrum of the AGN 
library. Further radiation components that may contribute to the 
IRAC fiuxes are ignored. The IRAC magnitudes of the library are 
computed using filter curves and fiux conversion as described in 
the IRAC Data Handbook. We find that most sources of the AGN 
library would be classified as AGN when the empirical mid IR 
selection criteria by [Stern et al.| ( |2005| ) is applied. Although there 
are some additional sources in the area near [5.8] - [8.0] >2, these 
are edge-on sources with high optical depth. 


presented in [Efstathiou et al.[ ( [20I3[ ) as in both models the mid 
IR emission is predicted to come from the conical region that 
is associated with the narrow-line region. The dominance of the 
mid IR emission of this conical region over the emission from 
the main body of the disk when the system is viewed approxi¬ 
mately edge-on is illustrate d in the image (to p middle) plotted 
in Fig.[^ We also note that [Xie et al.[ ( [2014] ) recently reported 
three other local ULIRGS with Spitzer spectroscopy which show 


similar SEDs to IRAS F10214-r4724. 


5. Conclusion 

We assume that dust in the AGN torus is distributed in a clumpy 
medium or in a homogeneous disk or as a combination of the 
two (i.e. a 2-phase medium). We have computed with a self- 
consistent three dimensional radiative transfer code the SED of 
such an AGN structure over a wide range of their basic parame¬ 
ters: the viewing angle, the inner radius, the volume filling fac¬ 
tor and optical depth of the clouds, and the optical depth of the 
disk midplane. We find that details of the cloud properties, such 
as cloud size, geometrical shape, or their density structure, as 
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Table 3. Fitting parameters of pure AGN models for representative objects. 



Model Parameter 

Derived Quantities 

Object 

Type 

Radius 

nn 

(lO^^cm) 

Eilling 
factor T] 

Cloud 

Disk 

^ V,mid 

Inc. 

6 

n 

Conv.^^ 

6 

Radius 

in 

(pc) 

Intrinsic AGN 
luminosity Lagn 
log(Lo) 

NGC 1365 

I 

7.8!^ 


14:51 

470+530 

^^'^-430 

22 22 + 10 
■^■^-14 


n 1 + 0.09 
'^•^- 0.06 

10 23+005 

+ +'•^■+-0.04 

NGC4151 

I 

4.8!^ 

22+75 

^-2 

40^39 

230+770 

-190 

43+® 

^■+-24 

2.08:0:15 

0 1+0-11 
+'• + -0.02 

10 35+^-^^ 
^'^•^^-0.02 

NGC 3081 

II 

3-4!f 


45+^ 

^^-44 

7+990 

-6 

8o:®7 

0.72:0;“ 

^•^- 0.02 

11 04+0-05 
+ +•++ -0.05 

NGC 5643 

II 


48!^ 

14+31 
^ -10 

400+600 

^+'+'-400 

67:11 

0.84:0;“ 

0 3+015 
+'•■+-0.2 

10 65+0-11 

iu . uj _ oo 2 

3C249.1 

I 

7 7^^ 
/./_2 



340:5“ 

6o:’ 

o-82:S:g 

4 8+51 

^•"-1.3 

13 57+^-^^ 

/-o 17 

3C351 

I 

7-8!^:^ 

14:“ 

i:i 

840:*“ 

43+17 

^■^-24 

9 + 0.04 

^- 0.05 

4 7 + 1-3 

^./- 2_7 

13 54+0-02 

3C324 

II 


9+74 

^-1 

4o:io 

290+700 

-170 

67:® 

- 0.14 

20 3+^-^ 

- 0.08 

3C356 

II 

10+5 

7+70 

-5 

4::° 

260+740 

^'++'-260 

6 T_f 

0 66+01® 

Q 9 + 4.6 
- 6.4 

1 9 q + 0.09 
^■^•^- 0.14 

F 10214^) 

+4724 

HE 

15:0-4 

4o:^^ 

6:^ 

910+90 

-600 

6 TJ , 

0 66+0-11 

- 

- 


Notes. Anisotropic luminosity conversion factor (Eq§. Inner radius of the dust torus /?in = V^agn/ 10^^- We do not 

estimate Ri^ and Lagn (Eq. for the lensed galaxy. 


well as the presence of clumps such as those detected by X- 
ray eclipse observations, have a minor impact on the AGN dust 
emission spectrum. 

We visualize the applied three-dimensional dust density 
structure of the AGN. The AGN torus emits anisotropically and 
generally the SED depends on the viewing direction. This hap¬ 
pens for wavelengths below the far IR peak flux. At longer wave¬ 
lengths the AGN emission becomes isotropic. We also show this 
in AGN images, where the scattering light (~ 1 pm) and the 
warm dust emission (~ 10 pm) is anisotropic, while the cold 
dust emission (>100yum) appears isotropic. The precise wave¬ 
length where the AGN becomes isotropic depends on the dust 
cross section. We And that type I sources emit more IR radiation 
than type IIs. When the disk material is included there is addi¬ 
tional hot dust in the system, which may increase the near IR 
flux for face-on views by an order of magnitude. 

AGN torus models often consider dust similar to what is as¬ 
sumed for the diffuse ISM. Because of the much higher density 
and stronger radiation environment in the AGN we favour fluffy 
grains. Their absorption cross section is generally larger than 
for ISM dust especially in the far IR and submillimeter being 
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a factor 10 higher at 1 mm. We show that for unresolved ob¬ 
servations, photon scattering within the beam may increase the 
detected flux and alter the wavelength dependence of the extinc¬ 
tion curve. This implies that there is no direct one-to-one link 
between the observed extinction curve and the wavelength de¬ 
pendence of the dust cross sections. Claims of detecting grain 
growth in AGN tori that are based on extinction measurements 
should therefore be treated with caution. 

The influence of the model parameters on the SED and their 
impact on the strength of the lOyum silicate band is discussed. 
The AGN library accounts well for the observed scatter of the 
10pm feature strengths and peak emission. The peak is centered 
near 10.4yum for ISM dust, whereas it is shifted to even longer 
wavelengths (~ 11.5 pm) for fluffy grains, and this is more con¬ 
sistent with observations. We cannot identify a striking need to 
postulate clumpy AGN torus models when explaining observed 
properties of the silicate band. However, arbitrary clump distri¬ 
butions of the AGN torus ease accounting for silicate features 
observed at various strengths and viewing angles. 

In the UV, optical and near IR, it may be necessary to add 
to the library SED a low luminosity component. This may be 
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Table 4. Fitting parameters of a combination of AGN and starburst models for representative objects. 



AGN parameter 

Starburst parameter 

Derived Quantities 

Object 

^in 

(10‘Vm) 

1 

^V,cl 

^ V,mid 

e 

o 

7?sb 

(kpc) 

T-sb 

log(Le) 

Ay 

(mag.) 

lOE 

(%) 

e 

(%) 

^SB 

^obs^' 

(%) 

Rin 

(pc) 

T-agn 

log(Le) 

3C249.1 

Q 7+0.2 

9+21 

^-0.5 

9+0.4 

^-0.3 

240+5‘« 

63!^ 

2.0![;® 

10.75!0;| 

126!^« 

83!’, 

0.82!«;« 

3 4+1-2 

3.9!j;’ 

13.28!«:“ 

3C351 

6.4!0;f 

12+11 

^^-10 

9+0.4 

^-0.4 

210+™ 

64!?o 

1 0+in 

^•^-0.3 

11.69tl 

00+61 

52!?2 

1 c 1 +0.04 
^•^^-0.05 

5.5!|’ 

4 9+2.1 
^•^-1.1 

13.75t«| 

3C324 

12.r0;« 

2+2 

^-1 

36!?6'8 

314+233 

68+^ 

1 7+1-2 

^-1.4 

11.88!0;3 

67!!^ 

46!r 

0 62+11-112 
-0.19 

13.8!“ 

Q 5+4.4 

13 76+11-111 

iO. /'J_o96 

3C356 

q 7+0.9 
•^•^-0.2 

1+0.8 

^-0.7 

9+2.2 

^-1.9 


70!’ 

1 5+1-2 

11.88!°;^ 

07 + 107 
^ ' -1% 

50!t° 

ri cq+0.09 
^•^"-0.06 

00 0 + 11.5 
00.0-17 7 

1 7+1-2 

^ • ' -0.4 

13 45+0-06 

i0.H-0_QQ3 

F10214 

-^4724 

3 6+°-^ 
•^•'^-0.6 

40+18 


34!f2 

27!'' 

1 7+0-3 
-1.0 

11.88!0;| 

92+51 

^^-39 

45!^ 

1 26+11-112 
^•^^-0.03 

19 1+2-^ 

- 

- 


Notes. Notation as of Table Ratio of the starburst and the apparent AGN luminosity. We consider for the starburst activity the library by 
I Siebenmorgen & Krugel| { [2007| ). 
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Fig. 14. Spitzer/IRAC mid IR color diagram. Observed A GN populate 
the ar ea in between the black lines at [3.6] - [4.5] > 0.3 ([Stern et al. 
|2005j ). Colours derived from the AGN library populate this area, lype 
I sources with viewing angle < 60® are shown as red circles and over¬ 
plotted are type II sources (> 60®) as blue circles . 


for type I sources beamed synchrotron radiation or for type II 
sources stellar light from the galactic disk. These photons es¬ 
cape the torus without interaction. Given a set of data points for a 
particular AGN, there is a simple procedure to select from the li¬ 
brary those elements which best match them. If the observations 
cover the full wavelength bands from the near IR to submm, one 
usually finds a library element that fits very well. This is demon¬ 
strated for a number of representative type I and type II objects, 
the famous hyperluminous infrared galaxy IRAS FI0214-1-4724, 
and 4 nearby Seyferts. In all cases, the parameters defining the 
element are in full accord with what is known about the object 
from other studies. Interestingly for the luminous AGN it is pos¬ 
sible to explain the SED of these objects with pure AGN models 


without a need to postulate starburst activity. In a second step 
we assume a combination of AGN and starburst activity with 
the goal of estimating an upper limit of the starburst luminosity. 
We find that the SED of the 5 high luminosity objects studied in 
this paper, require a starburst luminosity of typically 15% (with 
a range between 3 - 33%) of that of the AGN. In our model the 
starburst contribution to the far IR emission of the AGN can well 
be marginal. We also confirm that objects in the library span the 
range of mid IR colors as empirically established from observa¬ 
tions. 

Our library of AGN models promises to be very useful for 
interpreting the detailed ^o^i-Spitzer and post-Herschel SED of 
large samples of AGN and (ultra-) luminous infrared galaxies in 
general as well as high spatial resolution imaging observations 
with large ground-based telescopes and interferometers. The li¬ 
brary allows one to constrain the most fundamental properties of 
the dust torus which is powered by the central AGN. The SED li¬ 
brary also provides a utility that allows converting the observed 
(’apparent’) AGN luminosity, which is usually what is quoted 
in the literature, to the intrinsic (’true’) luminosity of the AGN 
source. The SEDs are accessible in a public library 
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